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Abstract 

We give a strong approximation of Rosenblatt process via transport processes and 
we give the rate of convergence. 
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1 Introduction 

Self-similar stochastic processes are of practical interest in various applications, including 
econometrics, internet trafhc, and hydrology. These are processes X = {X (t) : t > 0) whose 
dependence on the time parameter t is self-similar, in the sense that there exists a (self- 
similarity) parameter H € (0, 1) such that for any constant c > 0, {X (ci) : t > 0) and 
{c^X (i) : i > O) have the same finite dimensional distributions. These processes are often 
endowed with other distinctive properties. 

'Dedicated to the memory of Constantin Tudor 

^The author is supported by Dipuv grant 26/2009, ECOS/CONICYT C10E03 2010, MATHAMSUD 
09/05/SAMP Stochastic Analysis and Mathematical Physics Research Network. 

"''Partially supported by the ANR grant "Masterie" BLAN 012103. Associate member of the team Samos, 
Universite de Pantheon-Sorbonne Paris 1 



1 



The fractional Brownian motion (fBm) is the usual candidate to model phenomena 
in which the self-similarity property can be observed from the empirical data. This fBm 
is the continuous centered Gaussian process with covariance function given by 



2H 



t-sf). 



(1) 



The parameter H characterizes all the important properties of the process. In addition, to 
being self-similar with parameter H, which is evident from the covariance function, fBm 
has correlated increments: in fact, from ([1]) we get, as n ^ oo. 



when H < 1/2, the increments are negatively correlated and the correlation decays more 
slowly than quadratically; when H > 1/2, the increments are positively correlated and the 
correlation decays so slowly that they are not summable, a situation which is commonly 
known as the long memory property. The covariance structure also implies 



this property shows that the increments of fBm are stationary and self-similar; its immedi- 
ate consequence for higher moments can be used, via the so-called Kolmogorov continuity 
criterion, to imply that has paths which are almost-surely {H — e)-Holder-continuous 
for any e > 0. 

It turns out that fBm is the only continuous Gaussian process which is self-similar 
with stationary increments. This constitutes an alternative definition of the process. How- 
ever, there are other stochastic processes which, except for the Gaussian character, share 
all the other properties above for H > 1/2 (i.e. ([T]) which implies ([2]), the long-memory 
property, ([3]), and in many cases the Holder-continuity). In some models the Gaussian 
assumption may be implausible and in this case one needs to use a different self-similar 
process with stationary increments to model the phenomenon. Natural candidates are the 
Hermite processes: these non-Gaussian stochastic processes appear as limits in the so-called 
Non-Central Limit Theorem (see [2], [5], [K]) and do indeed have all the properties listed 
above. While fBm can be expressed as a Wiener integral with respect to the standard 
Wiener process, i.e. the integral of a deterministic kernel w.r.t. a standard Brownian 
motion, the Hermite process of order g > 2 is a gth iterated integral of a deterministic 
function with q variables with respect to a standard Brownian motion. When q = 2, the 
Hermite process is called the Rosenblatt process. This stochastic process typically appears 
as a limiting model in various applications such as unit the root testing problem (see \20\ ) 
or semiparametric approach to hypothesis test (see [IQJ). On the other hand, since it is 
non-Gaussian and self-similar with stationary increments, the Rosenblatt process can also 
be an input in models where self-similarity is observed in empirical data which appears to 
be non-Gaussian. The need of non-Gaussian self-similar processes in practice (for example 
in hydrology) is mentioned in the paper [17] based on the study of stochastic modeling for 



E [{B^ (n) - B^ (1)) B^ (1)] =H{2H- 1) n^^'^ + o {n 



2H-2 



(2) 




(3) 
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river-flow time series in [TT]. Recent interest in the Rosenblatt and other Hermite processes, 
due in part to their non-Gaussian character, and in part for their independent mathematical 
value, is evidenced by the following references: [T], [H], [Tl], [TS], [TU] . 

In this paper we will give a strong approximation result for the Rosenblatt process 
by means of transport processes. It is also interesting from the theoretical point of view 
since all the approximation results for the Rosenblatt process known in the literature are 
in the weak sense ([5], [IB]). 

Our work is a natural extension of the strong approximation results for the Brownian 
motion and for the fractional Brownian motion. The study of the convergence of transport 
processes to the Brownian motion has a long history. We mention the works ([1], [9], [8]) 
among others. More recently, due to the development of the stochastic analysis for fractional 
Brownian motion, the need of simulating the paths of this process led to the study of the 
strong approximation of the fBm by means. We refer to [7] for such an approximation in 
terms of transport processes and to [6] or [15j for related works. 

Our paper is organized as follows. In Section 2 we give some preliminares on multiple 
integrals and Malliavin derivatives. In section 3 we describe the approximating processes 
and prove the convergence to Rosenblatt process. 

2 Multiple Wiener-Ito Integrals and Malliavin Derivatives 

We start by introducing the elements from stochastic analysis that we will need in the 
paper. Consider Ti a real separable Hilbert space and {B{ip), (p € %) an isonormal Gaussian 
process on a probability space (0,^, P), which is a centered Gaussian family of random 
variables such that '&{B{ip)B{ilj)) = {ip,ip)-}{. Denote by I„ the multiple stochastic integral 
with respect to B (see ^13j). This /„ is actually an isometry between the Hilbert space 
^®"(symmetric tensor product) equipped with the scaled norm -^r^i II ' ll'H®" ^^'^ Wiener 
chaos of order n which is defined as the closed linear span of the random variables Hn{B{ip)) 
where if G 7^, II'^II'h = 1 and Hn is the Hermite polynomial of degree n > 1 




The isometry of multiple integrals can be written as: for m, n positive integers 



E(I„(/)/^(<7)) 

W.{InU)Im{g)) 



n\{f,g)'^»n if m = n, 
if m 7^ n. 



(4) 



It also holds that 
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We recall the following hypercontractivity property for the norm of a multiple stochastic 
integral (see [12, Theorem 4.1]) 

E|/„(/)|2-<c™(EI^(/)2)"^ (5) 
where Cm is an explicit positive constant and / G 'H®'". 

In this paper we will use multiple stochastic integrals with respect to the Brownian motion 
(Bm) on M as introduced above. Note that the Brownian motion on the real line is an 
isonormal process and its underlying Hilbert space is T-L = L^' 



For every ^ < H < 1 the Rosenblatt process (^i^)t(:[o,T] could be defined as follows, 

Xl' = ciH)h{gti-)) (6) 

where for every t £ [0, T] 



gt{yi,y2)= {u-yi)j^ (n - y2)+ du. (7) 

JyiVy2 

The constant d{H) is a normalizing constant which ensures that E(X/^)^ = t"^^ for every t G 
[0, T]. This constant can be explicitly computed but it has no interest for our investigation. 
It can be proved that the process X^ is self-similar with stationary increment and has the 
same covariance ([T]) as the fBm. Moreover it satisfies properties ([2]) and ([3j). 

3 Strong convergence to the Rosenblatt process 

The Rosenblatt process (-^/^)te[o,T] defined above can be written as an iterated double 
integral in the following way 

Xl^ = c{H)[ f ( f\s-xi)f'\s-X2)J'^ds)dB{xi)dB{x2), t e [0,T] (8) 

where i? is a Wiener process on the whole real line and the Hurst parameter H belongs to 
the interval 1). The process X^ is H self similar with stationary increments and it has 
the same covariance as the fractional Brownian motion. 

We will separate X^ into three terms. For every t E [0, T] 

r rO rO / rt 

X? = c{H) 



[ [ ( [ {s-xi)f-\s-X2)^-^ds\dB{xi)dB(x2) 
+ [ [ ( [ (s- xi)^-\s - X2)^~^ds) dB{xi)dB{x2) 

J-oo Jo \Jxi J 

+ [ [ ( [ {s-xi)"-\s-X2)^'^ds)dB{xi)dB{x2) 

Jo J-oo \Jx2 J 
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+ f f ([ {s-xi)^~\s-X2)^~^ds)dB{xi)dB{x2) 

Jo Jo \JxiVx2 J 

:= x/'^ + 2X,2'^ + xf'^, (9) 

note that the second and the third integrals are actually equal, for that reason the term 
X^'^ appears twice . We will treat separately the third terms above since they have different 
behavior which comes from the singularity of the integral appearing in their expression. 



3.1 Transport processes 

For each n = 1, 2, . . ., let (Z("^(t))t>o be a process such that Z^^\t) is the position on the 
real line at time t of a particle moving as follows. It starts from with constant velocity +n 
or —n, each with probability 1/2. It continues until a random time ri which is exponentially 
distributed with parameter n^, and at that time it switches from velocity in to and 
continues for an additional independent random time T2 — ri which is again exponentially 
distributed with parameter n?. At time T2 it changes velocity as before, and so on. This 
process is called a (uniform) transport process. Griego, Heath and Ruiz-Moncayo [9] showed 
that converges to Brownian motion strongly and uniformly on bounded time intervals, 
and a rate of convergence was derived by Gorostiza and Griego in [8] as follows. 

Theorem 1 There exist versions of the transport processes on the same probability 
space as a given Brownian motion {Bt)t>o such that for each q> 



P I sup \Bt - Z["^ \ > Cn~^/^{lognf/'^ ) = o{n-'^) as n ^ oo, 

. a<t<b 



where C is a positive constant depending on a, b and q. 
Let (Xf ) 

te[o,T] a Rosenblatt process. With a < fixed, we consider the following 
Bm's constructed from the Bm B in 

1. (-Bi(s))5gjo T] ' restriction of B to the interval [0,T]. 

2. (-B2(s))^<^<Q, the restriction of B to the interval [a, 0]. 



3. B^is) 



sbC-) ifse[i,o), 

if s = 0. 



Let us define now the transport processes that will intervene in our main results. By 
Theorem [H there are three transport processes 

iz'r\s))o<s<T, {Zt\s))a<s<0, and {zt\s)) (10) 
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such that for each q > 0, 

p( sup >C»n-i/2(logn)5/2 ) =o(n-'') as n ^ oo, (11) 

\bi<t<c, J 

where h, Ci, i = 1,2, 3, are the endpoints of the corresponding intervals, and C^*^ is a positive 
constant depending on 6j, Cj and q. 

3.2 Strong approximation 

We will approximate successively each summand X^'^, X^'^, X^'^ from ([9| in the strong 
sense by processes construct in terms of the transport processes Z;^"^ introduced 

above. Let us start with the summand X^'^ . Using Fubini theorem, we can express it as 

rO rO 



where 



X]'^ = c{H) ds / dB{xi)dB{x2){s - xi)--^{s - X2)~'^ 

J J ~oo J —oo 

= c{H) j^ds(^j {s-x)^-^dB{x)^ (12) 
= c{H) f dsiY^^'^f, tG[0,T] 

JO 

Y}^^ = {s-x)--^dB{x), se[o,r]. (13) 



Remark 1 Notice that integral J^^o^s — x) 2 ^dB{x) is well-defined in L'^{Q) as a Wiener 
integral for every s > since 



{s-x)^~^dB{x)\ = I {s-x)"-^dx = — 

\J—oo J J— 00 1 



_J_„2H-1 



The situation will be different when we treat the summand X^'^ . This is one of the reasons 
to decompose the Rosenblatt process into several parts. 

Let < max ^ \J^{^ , 2H+2 ) < P < 1/2 be fixed (note that < \ since H < 1 

and 2H+2 ^ h because H > denote in the sequel by 

£^ = n (14) 

and by 

a„ = n~*^2~'^)(logn)5. (15) 
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We will use the notation 



ll^lloo.M = sup \Ys\. 

a<s<b 

When the interval is of the form [0, T] we will use the shorter notation t] ■~ ll^lloo T • 

We will denoted by C a generic strictly positive constant that may depend on a, T, H,p and 
may change from line to line. 

Let us give a different expression for the process Y^'^ . 
Lemma 1 Let Y^'^ be the process defined by il3\) and a < fixed, then for every s G [0, T] 

= fs{a)B2{a)- [ d,fs(-)^B3{u)du- r d,fJ-)XMu)du 

r-Sn rO 
+ fs{x)dB2{x)+ iUx) - fs{x - en)]dB2{x) 

J a J —Sn 

+ [ fs{x - en)dB2{x) (16) 

■J ~£r,. 



where dxfs denotes the derivative of the function fs{x) = (s — x)^/^ ^,s > x with respect 
to its second variable (even when this second variable is not denoted by x). 

Proof: We can write, for every t € [0, T] 

y}'''= r fs{x)dB{x)+ f fs{x)dB{x). (17) 

We express the first Wiener integral above as an integral with respect to ds. Since by the 
Holder continuity of B, 

lim fs{b)B{b)=0, 

b—^—oo 

by integration by parts and putting x = 1/u, 

r fs{x)dB{x) = fs{a)B{a) - T dxfs{x)B{x)dx 

J — oo J ~oo 

= fs{a)B{a) - f dxfs (-) \b (-\ du 





fs{a)B2{a) - I djs (-) ^B^{u)du. (18) 



'1/a 

By pT|) and (fT8|) . for every s € [0, T\ we have the result. ■ 

We first approximate the process (^s^'^)se[o,r] (iu the strong sense (fTTI) ) by stochas- 
tic processes constructed from transport processes. Basically, in the expression of Y^'^ , we 
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replace the Brownian motions by their corresponding transport processes. The approximat- 
ing processes to Y^'^ is defined as 



l,H,n 



l/a 



+ f fs{x-en)dzl^\x), se [o,r]. 

J —e-n. 



(19) 



We state the result concerning the approximation of Y^'^ . Its proof follows the ideas 
of the proofs in [7] but the context is technically more complex. Note that the singularity 
of the integrand {s — x)^~^ at s = x does not allows to use directly the results in [7j and 
the arguments of the proofs must be adapted to fit in our context. 

Proposition 1 Let Y^'^ and Y^'^'^ be the processes defined by /il3\) and I119\) . respectively 
and let an given by il5\). Then for each q > and each (3 such that < < /3 < ^, 



P sup S^--^/^ \Yl,H _ y^l,/f,n| > \ ^(^-9) aS U ^ OO. 

\Q<s<T J 

Proof: From (jl6p and Lemma [1] we have 



(20) 



ft{a)B2{a) - ft{a)Z^^\a 



l\ 1 



l/a \UJ U'^ Jl/a \uj U-^ 





dxfs ( - ) -KB2,{u)du 



dxfs { - ] -Z'i'\u)du 
fs{x)dB2{x)- I fs{x)dZ^^\x) 

J a 

+ f fs{x-en)dB2{x)- C fs{x-£n)dZ^2\^) + [/. (x) - /. (x - e„)]dS2 (x) 
By Lemmas El [3l [H [5l [6] and [7] below we have the result. ■ 

Lemma 2 Let Z^^ be the process defined by II 10]) . Then for each q > there is C > such 
that 



Ii := P sup 

V 0<s<T 



fs{a)B2{a) - fs{a)Z^''\a) > Ca„ j = 0(7^-") as n^oo. (21) 
Proof: It holds, for fixed a < 0, 



fs{a)B2{a) - fs{a)Zt\a)) < ||i?2 - ^^"^ ||oo,[a,oi(^ " a)^/^-^ 
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< \\B2 



-'2 llcxD 



then (recall that C is a generic strictly positive constant that may depend on a, T, H) 
h < p(||B2-^?^||oo,M](-«)''/'-'>Ca„) 



< P[\\B2 - ^^"^||oo,[a,o] > Can) = oin''^). 



Remark 2 The conclusion of Lemma\^ is clearly true if we add the factor s^~~ after the 
supremum. This remark is also available for the following lemmas and we will not mention 
it at each time. 

Lemma 3 Let Z,"^ be the process defined by ^0\) . Then for each q > 0, 



h := P{ sup 5^-^/2 
\ 0<s<T 



'"dj, (-) l^Bs{u)du- [ '"d^fs (-) ^Zi''\u)du 

l/a \U/ U Ji/a \U/ U 



= o(n as n oo. 
Proof: Putting z = 1/u and w = s — z, 



l/a 



dxfs ( - ) -^B^{u)du 



< 



< 



< 



< 



\\B3 


- ^3 \ 


\\B3 


- ^3 1 


\\B3 




\\B3 


- ^3 1 


\\B3 


y(n) 1 
- ^3 1 


\\B3 


y{n) 1 
- ^3 1 


\\B3 


yin) 1 
- ^3 1 


\\B3 


yM 1 
- ^3 1 



oo,[l/a,0] 



'l/a 

,[l/a,0](l - H/2) 



l/a \UJ W 

dccfs (-) ^B:,{u) 



(u)du 
du 



l/a i-u)^ 



- l/n)^/2-2dn 



oo,[l/a,0](l - H/2) / {-Z){S - z f/'-'dz 

J-l/e-a 
i-s+l/en 

loo,[l/a,0](l - H/2) / {W - s){wf/'-'dw 

J s—a 

/s+l/e„ 
_ —a 

loo,[l/a,0]i^[(« + IK)^/^ - - a)^/^] 



H/2 

loo,[l/a,0] {T + \/enf^'^ 



Wi/aM-^2^/\T^I^ + {en)-^l% 



> Car 



(22) 



then 



<p(||i?3-^f^||oo,[l/a,0] >Ca 

+P (1IB3 - 4"^||oo,[i/a,o] > C{enf/^a 

< o{n-'^) + P (1IP3 - ^l"^lloo,[i/a,o] > Cn-V2+/^(i~^^)/(i-^/2)(io; 
= o(n-«). 



<P( ||i?3-Z(")||^,[l/, 



+P I IIP3 - Z3 - ||oo,[l/a,0] jjj^ 



n 



The following lemma explains one of the conditions imposed on (3 in the statement 
of Proposition [TJ Another restriction comes from Proposition 4 later. 

Lemma 4 Let (1 - H/2)/{3 - 2H) < (3 <\. Then for each q>0 



h-=P{ sup 

\ 0<s<T 



> a. 



o{n as n ^ 00. (23) 



Proof: For every s £ [0,T] and u € [— e„,0) we can write 



1 



(^/2-l). 



l/n)^/2 



-2 



l-H/2 (l-us 



-u 



-u 



H/2-2 



< {1- H/2){-u)-^-^/^ 



(24) 



By (j24p and the pathwise Holder continuity of the Bm P3 there exists a random 
variable Y (having all its moments finite) such that for any 7 < 1/2 — H/2, 



/1\ 1 

dxfs - —Bj,{u)du 



<Y I (1 - H/2){-uY^-^/^{-uf/^-^du 

1 ~ -^/2 \ll2-H/2-y 



Y 



1/2 -i?/2 -7 



By Chebyshev's inequality, for r > 0, 
h 



<p\y n 



1/2- H/2 --i 



„ l/2-g/2-7 
>■ l-il/2 



> Or 
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= P(CY> n^ilognf/') < 

where k = -(1/2 - /3) + /3(l/2 - H/2 - 7)/(l - H/2). Taking 

(1 - H/2)/{3 - 2H) < (1 - F/2)/(3 - 2// - 7) < /3 < 1/2, 
then K > 0. For g > there is r > such that q < rK, then 

hm n^/s = 0. 



Lemma 5 Let he the process defined by i TfOj) . Then for each q > 0, 



I4 := PI sup 

\0<s<T 



fs{x)dZ^''\x) 



> Can 



(25) 



(1 - H/2){s - xf^^-^B2ix)da 



fs{x)dB2{x 

= o{n~'^) as n —7- 00 
Proof: By integration by parts, 

/ fs{x)dB2ix) = fs{-en)B2{-en) - fs{a)B2{a) 

J a 

and 

r" f,{x)dZ^i'\x) = M-en)zt\-en) - fsia)zi''\a) 

J a 

-/ (l-^/2)(s-x)^/2-2z^")(x)d; 
J a 

then 

/ fs{x)dB2{x)- I fs{x)dz!^"-\x 

J a J a 

< \\B2 - 4'^||oo,[a,o] |/.(-en) + fs{a) + J "'(1 " H/2){s - xf/^-^dx 
<\\B2- |U,[„,o] {{s + + {s- a)^/2-i + {s + ey-' - {s - a)^''-'} 

= ||i?2-^^"^||oo,[a,0]2(s + 

< 11^2 - z(")|U,[,,o]2(en)''/'"' = \\B2 - z(")|U,[,,o]2n^ 
Consequently, 

h <P [\\B2 - ^ ||oo,[a,0]2n^ > Can) 

= P (\\B2 - Z?^||oo,[a,o] > Cn-'/Hlognf) = o(n-'') 
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Lemma 6 Let z''^^ he the process defined by [W\) . Then for each q > 0, 



sup / fs{x - en)dB2{x) - / 

\Q<S<T J-En J-t 



:= P 

= o{n^'^) as n ^ oo. 



> Car. 



(26) 



Proof: By integration by parts as before and taking into account that -62(0) = (0) = 
we can express the two integrals in the statement as 

/ fsix- en)dB2{x) = -fs{-2en)B2{-en) - (1 - H/2){s + e„ - x)""'^'^ B2{x)dx 



and 



f,{x-en)dZ^i'\x) 



- / (1 - H/2){s + £n - x)^/^"^Z^-^\x)dx. 
J -en 



then 



f fs{x-en)dB2{x)- r fs{x - en)dZ^'"\x) 

C ||i?2-Z^")||^,[„,0]{/s(-2en) + f (l-i?/2)(s + 5„-x)^/2-2dx| 

C \\B2 - z(")||^,[„,o] {{s + 2£„)^/^-i + (. + Enf'-^ -{s + 2e„)^/2-i} 



\B2 - Z 



2 l!oo,[a,0] 



Consequently, 



< P 



(||i?2-^^"^||oo,[a,0p''>Can) 



P I ||i?2 — Z2 Hoc 



[,,0] > Cn-i/2(log 



n 



o[n 



Finally, we prove our last auxiliary approximation result. Here we need to add the 
factor s^~~ which appears in Proposition [TJ This is due to the singularity of the derivative 
of fs{x) with respect to x. s^~^ . 
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Lemma 7 Let (1 - H/2)/{3 - 2H) < /3 < i Then for each q>0 



I -1 

P sup S 2 
\ o<s<r 



r [/s(x)-/,(:E-e„)]di?2(x: 

J —En 



> Car, 



oin 



as n ^ oo. 



(27) 



Proof: We first write the difference fs{x) — fs{x — En) as an integral and we use Fubini 
theorem. We obtain 

l-O r-S+En—X 

[fs{x) - fs{x - en)]dB2{x) = / / {l-H/2)u''/^~^dudB2{x) 

-Sn ^ Sn S — X 

/ dB2{x)du+ / 

J S~U J S+En 

S + En 

u"/^-^[B2iO) - B2is - u)]du 



(1 - H/2) 



dB2{x)du 



(1 - H/2) 

s+2e 



rs+ZE„ 

+ / u^/^-^[B2{s + en-u)- B2{-en)]du 

J S+En 



The Holder continuity of the Wiener process B2 implies for every < 7 < i 



[fs{x) - fs{x - en)]dB2{x] 



< (1 - H/2)Y 

< (1 - H/2)Y 

< (1 - H/2)Y 

and consequently 



i+En /•S + 2En 



S+En 



+£71 



S + 2Er. 



I -I r. 

P sup S 

V 0<s<T 



- fsix - en)\dB2{x) 



> Can ] <P CYel/^^^ > 



a. 



and the result follows by analogous arguments as in proof of Lemma [H 
Remark 3 In particular Proposition [7] implies that 

y " o<s<r 
by using Borel-Cantelli lemma. 



sup S^-'^W^^H _yl,H,n\ 
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We finish the strong approximation of the term X^'^ appearing in the decomposition 
of the Rosenblatt process in dSD. By we define for every n and Y ' ''^ given by 

m 

XhH,n ^ f ^Y},H,n^2^^_ (28) 

Jo 

We have the following. 



Proposition 2 Let X^'^ be given by flgj) and /3 G y \-^h ' ^) /^2;e(i. Define X^'^'^ by 
). Then for any 7 such that < 7 < /3 and f3 + j < 1/2, 



Proof: Using the fact that — B'^ = {A — B)'^ + 2B{A — B) we can write, for every 

t G [o,r] 

Xl^^'^^-Xl^^ = c{H) [\{Y}'^'-f - {Y}'''f)ds 

Jo 

= c{H) f [(y/'^'« - Y}^^f + 2Y}^^{Y}'^^- - Y}^^)] ds 
Jo 

and hence 



sup < c{H) sup [ \{Y}'^'--Y}'^f + 2Y,''^{Y}'^'--Y}'^)\ds 

te[o,T] te[o,T] Jo 

= c{H) / \{Y}'"^'^ -Y^^^f + 2Y}^^{Y}'"^'^ -Y}^^)\ds 
Jo 

(■T rT 

< c{H) / (Yy''' -Y}'")'^ds + 2ciH) / IY^^^Wy}'"^"" -Y} 

< c sup s2-^(y,i'^'" -y,i.^)2 

se[o,T] 

T 

+ 2C [ sup s^-f |y/'^'"-y,i'^|. 

^0 se[o,T] 

We used above the trivial inequality > Ca„) < > Ca„) for any random 

variable X. We will get (C denoted a generic strictly positive constant depending on T, 
that may change from line to line) by Proposition [1] 

P ( - Xl'^lloo.T > Cn-(V2-/3-7)(iog„)5/2 



< P sup s2-^|y,i'^'"-y,i'^|2>Cn-(i/2-/?-^)(logn)5/2 
\se[o,T] . 
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pf sup ^i-f |y,i'^'"-y,i'^| rsf-i|y,i'^|d5>Cn-(i/2-/3-7)(iogn)^ 
\se[o,T] Jo J 

o{n-'i) + pl r s^-^Y}'"\ds sup s^-f ly^i'^'^^-y,!'^! > C7n'(i/2-/3-7)(iogn)5/2 I 
\Jo se[o,T] / 



(29) 



A= I s^-^\Y}'^\ds and r= sup s^"^ ly.^'^ - y.^'^'"! 



We apply Lemma [5] below with 
rT 

lO 0<s<T 

We note first that 

'''^ s^-'^\Y}'^\ds < c{H) r 
'0 Jo 



f'l f 1 I 

E s^-'^\Y}'^\ds < c{H) ds s^-^ {EiY^'^fy ds 
Jo Jo 

< c{H) [ ds s^-^ ( [ (s - x)^~'^dx 





= c{H) [ ds s^~^s^ = c{H)T"-^ 
Jo 

and thus the random variable A is almost surely finite. We obtain by (|29p and Remark [3] 



< P lim {A sup s^'- |y/'^ -y,i'-^'"| > Cn-(i/2-/3-7)(iog„)5/2| 



< P I lim { sup s^-- \Y}'^ - y/'^'"| > Can} = 0. 



The following lemma has been used in the proof of Proposition [2l 

Lemma 8 Let A and T be random variables with A an almost surely finite. Then for every 
7 > 



P ( lim {AT > Cn-(i/2-/3-7)(iog^)5/2|\ 
< P f Ik^^ {P > C7n-(^/2-/3) (iQg ^yo/ 



with C a generic strictly positive constant. 
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Proof: We prove the following inclusion 

(at > n-(i/2-/3-7)(iogn)5/2l c Ik^I |r > n-(i/2-/3)QQg^)5/2| _ 

Since 

ue Ih^ f^r > n-(i/2-/3-7)(iog^)5/2l ^ jjj^ \{An-^)r > n-(^/2-/3)(iQg^)5/2l 

= n- 1 ur=„ {(ylA;-^)r > A:-(i/2-/3)(iogfc)5/2| , 
then for all n > 1, 

CO G U^=„{(AA:-T)r > A;-(V2-/3)(iogfc)5/2| ^ 

and since ^ is an almost surely finite random variable, there is iV = N{uj) such that for 
ah n> N, An~^ < 1, then u £ {r > ^"(^/^^^^(log A;)^/^} and the conclusion follows 
easily. ■ 

Let us handle now the term X^'^ appearing in Q. We will decompose it as follows: 

t ft ft 



X^'" = c{H) [ I dB{xi)dB{x2) I ds{s - xi)^'^{s - X2)^'^ 

Jo Jo J(xiVx2+en)/\t 

ft rt p{xi\JX2+en)/\t ^ ^ 

+c{H) / / dB{xi)dB{x2) / ds{s - xip-^{s - X2p~^ 

Jo Jo JxiVx2 

:= c{H)iF,- + G^) (30) 



where e„ is given by (jl4p . 

Remark 4 For the term X^'^ we cannot use Fubini theorem (as in the case of X^'^ ) 
because fQ{s — x)~~^dB{x) is not defined as a Wiener integral in L^(J7) since the function 
(s — is not integrable on [0,s] with respect to dx. 

For every t € [0, T] the summand can be written as 

Fj" = / / dB{xi)dB{x2) [ ds{s-xi)^'^{s-X2)^'^ 

Jo Jo J {Xi_\JX2+Sn)/\t 

/ dB{xi)dB{x2) ds(s + - xi)~"^(s + e„ - X2)~"^ 

Jo J(xiVx2) 
ft-En / fs ^ \ 2 rt-en 

j U is + en-x)--UB{x)] ds = J {Y.^'^fds (31) 



where we denoted by 



rs 

Y^'^ = / (s + e„ - x)f-'^dB{x) for < s < T. (32) 
Jo 
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Note that the process Y^'^ depends on n. But we prefer to use the notation Y^'^ without 
n in order to keep the coherence with the other terms treated before and in the sequel. 

Let (-Bi(s))sg[o,T] the restriction of the Wiener process B to the interval [0, T] and 

let Z^""^ be the corresponding transport process defined in (jlOp that converges to B in the 
strong sense (fTTIl . Then, 



Ys^'" = [\s + en-x)^-^dBi{x) 
Jo 



(33) 



and we define 



y/'^'" = / {s + £n- x)T~^dz["'\x) for < s < t. 

^0 

We will show first that y3,H,n jg strong approximation of Y^'^ . 

Proposition 3 Let Y^'^ , y3,H,n given by (3^, [34\ ) and l[T5\) respectively, 

P I sup |y/'^'" - y/'-^l > Ca„ I = o{n~'i) 
V^e[o,r] J 

for each q > and for (3 G (0, i). 

Proof: After integrating by parts, we can write, for every s G [0, T], 

B,{s)-Z[^\s) 



(34) 



\J-s I ^ 



f'S 

+{l-H/2) / (s + Sn-x)^-^ Bi{x)-Z[''\. 
Jo 

< 2ef'^\\Bi - ^||oo,T = 2n^\\Bi - zi")||oo,T 



dx 



using the choice of e„ and hence by (fTTj) . 



3,//,n _ y3,H 



p iiy 



for every g > by pT|) . 



,,T > Cq„) < P ||Pi - Z| 



|oo,T 



> C7n-i/2(logn)5/2^ =0(71"'') 



We will introduce now the approximation processes that will converge to X^'^ . Let 
us denote, for every t G [0,T], by 



X 



3,H,n 
t 



c{H) fiY^^^'-fds. 
Jo 



(35) 



The part X^'^ is approximated as follows. 
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Proposition 4 For < max \ ^_2h -> 2+2H j < < ^ fixed, let X^'^ and X^'-^'" defined 
by Ii30\) and i35\) respectively. Then for every 7 such that < 7 < /3 and 7 + /3 < ^, 

Proof: By (f30l) and (f3T]). 

Jo 

= c{H) f [(y/'-f^ - y/.-f^'")2 + 2y/'^(y/'^'" - y^^'^)] ds - c{h)gi 

Jo 

we have the bound 

Jo 

(■T 

+ 2c{H) / \Yf^^\\Yf^^ -Yf^^^''\ds + c{H) sup 
JO te[o,T] 

The first two summand in the right hand side above can be treated as in the case 
of X^'^ in Proposition [2j We note, in order to apply Lemma [HI we need to notice that 



E 



j^Ys' ds < C with C not depending on n (Lemma [8] can still be used although the 
process y^.-ff depends on n). Let us handle the term G". We will actually show that 

Vp( sup > Cn-(i/2-/3-7)(iog^)5/2) ^ ^_ (3g) 

te[o,T] 

which will imply that 

P fc{|jG"||oo,T > Cn-(l/2-/3-7)(iog^)5/2|\ ^ ^ 



For every t G [0,T] we have, using the fact that the integrand is symmetric with respect to 
the variables xi and X2, 



= 2 / dB{xi) / dB{x2) / ds{s-xi)"^is-X2)^~^ 

Jo Jo Jxi 

(>(t-e„)VO i-xi rxi+£n ^ ^ 

= 2 dB{xi) / dB{x2) / ds{s - xi)^-^{s - X2)^~^ 

Jo Jo Jxi 



--It 



J{t-e„)yo Jo Jxi 

^l,n^^2,n_ 
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Note that the mapping 

xi ^ dB{x2) / ds{s — xi)^~'^{s — X2)^^'^ 

Jo J XI 

is adapted with respect to (the filtration generated by the Wiener process B). Then the 
process {G\'"')t is a martingale for every n. Then we have, taking = n~(^/^~^~'''^(logn)^/^ 
and using Doob's inequality, 



P\ sup GJ'" > Cn-(i/2-/3-7)(iog„)5/2 <a-PE sup 



^te[o,T] 



G 



l,n 



Xl+En 



(iS(xi) / dB(j;2) / ds{s - xi)^~^{s - X2)^'^ 

^0 Jxi 



with C allowed to depend also on p in this proof. Note that the random variable 
j|(r-e„)vo ^jx^ ^^(^2) J^l^'"-" ds{s - xi)^-^{s - X2)^-^) is a multiple integral of 

order two. Therefore, by the hypercontractivity property ([5]), it is not difficult to see that 



P sup G,^'">Cn-(i/2-/3-7)(iogn)5/2 

\t6[0,T] 



< Ca-P 



dxi I dx2 ( / 

^0 \Jxi 



ds{s - Xl) ^ ^{s — X2)^ ^ 



Let us first compute the integral with respect to ds. By making the change of variables 



s—xi 

S — X2 



with ds 



X\—X2 



dz we get 



f 

Jo 

f 

Jo 



rT-£n j-Xl / fXi_+£n 

dxi I dx2 I / ds{ 
Jo \Jxi 



s — xi) 2 -^(s — X2) ^ ^ 



rT-En rxi 

dxi / dx2{xi — X2)'^^^'^ 

Jo 



We separate the integral dxidx2 into two regions: when X1 — X2 < £n and when X1 — X2 > £n- 
The above term will be bounded by 



dxi 



fXl 



'{xi 



dx 



,)vo 



2{Xl - X2 



,2H-2 



H_ 
Z 2 



-1^ 



dz 
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r dxi r dX2{xi - X2f"-^ { f Z^-^{1 - z)-^d^ 

Jo J(xi-en)VO \Jo J 

pT /•(j;i-£„)VO / r „ \ 

1^ dx,£ dx2(xi-x2r-2 n^""""^%f-\i-z)-^ciz 



+ 

< c{H) I dxi 



XI 



dx2{xi - X2) 

Xl-Cn 

(xi-e„)VO 



2H-2 



+c{H) f dxi [ ' ' dx2{xi-X2f"-^F^{H/2,H,H/2 + l,l] 
Jo Jo 

where F1{H/2,H,H/2 + 1, 1) is the incomplete beta function and hence 



En + Xi- X2 



H 



P { sup Gj^'" > ca„ ) < Ca-Pe^n" 
V*e[o,T] ) 



and the series 5^„£n^^^ oirJ' is finite if 

Note that (5 > imphes that /3^^ - - P — j > iov small 7 > 0. By choosing p 

large enough we obtain that 



V P ( sup Gj^'" > Can ] < 
\telo,T\ J 



00 



for every /3 € (0, |). 

Let us handle now the term denoted by G^'". We have 



P I sup Gf''' > Can I < oc-P^ sup 



,te[o,T] 



*6[0,T] 



G 



•2,n 



G 



we will use Garsia's lemma. To this end we need to 



In order to control Esupjgjo,T] 

estimate the LP norm of the increment Gj " — Gg'" when t is close to s. Note first that, by 
the change of variables z = we have 

° S—X2 

Gj'" = / dB{xi) r dB{x2)\xi - X2\"-^ z^-\l - z)-"dz 

Jt-£n Jo Jo 
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and for t,s € [0, T] such that s > t — e„, by the isometry of multiple stochastic integrals (jH) 

2 



E 



r-t r-xi 

c{H) / dxi / (ix2[a;i— X2 



- z)~"dz 







t-e JO 



XI 



-2c{H) I dxi I dx2\x, - X2|2^-2 I / z^-Ul - zr"dz 1 I / zf -1(1 - z)-"dz 



and by majorizing Jq* z 2 ^(l — ^) ^dzhyj^z'^ ^{1 — z) ^ dz = A — H) we oht 
the bound 



am 



E 



^2,n _ (^2,n 



2 /"^l 

< c(//) / dxi / (ix2|xi— rc2 

it-e„ Jo 



2H-2 



l^t pxi rt pxi 

+c{H) I dxi / dx2\xi - X2\'^^^'^ + c{H) / dxi / dx2\xi - X2\^^~'^e 

J S—Cn Jo J S—En JO 



The hypercontractivity property of multiple integrals dS]) implies that for all t, s 



E 



/-~(2,n ^2,\ 



<c{p,H){E 



/~(2,n ^2,\ 



2\ 2 p 

< Cel. 



where C is a constant that depend on H and p. Finally, by Garsia lemma (see e.g. |13j , 
Appendix A. 3) for p > 2 

E sup < Cel (37) 



0<t<T 



for every 7 such that < 7 < | — 1. The bound (f37|) implies, using Markov's inequality 
and taking suitable p large enough, that 



P sup G^'" > cdn I < 00 



(38) 



1-^ 



due to the fact that /3 > 3^2/7 ^'^cl this finishes the proof. 



Let us finally treat the summand X"^'^ in ([9]). Its approximation will be a mixture 
of the approximations of X^'^ and X^'^ . We have 



/ dS(xi) / dS(x2) / {s - xi)^-'^{s - X2)^~'^ds 

J —00 Jo J X2 
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/ dB{x2) / Y}'^{s-X2)-~^ds 

Jo Jx9. 



with Y^'^ given by (|13p . To avoid the singularity of the integral with respect to ds at 
s = 0:2 we will decompose this integral into two parts. In this way we can write, for e„ the 
sequence converging to as n — >■ 00 chosen before 



Xt'^ = [ dB{x2) [ Y}'"{s-X2)^-^ds 

Jo J{x2+en)At 
rt r{x2+en)/\t 
+ / dB{x2) / Y}^''{s-X2)"^ds 

Jo Jx2 

ft fS~e„ ^ ft f(x2+en)/\t ^ 

/ dsY}'" {s-X2)^-^dB{x2)+ dB{x2) Y}'^ {s - X2)^-^ds 

Je„ Jo J(t-en)VO Jx2 



^(t-e„)VO f{x2+en)At ^ 

+ dB{x2) Y}'^{s-X2)--^ds 

Jo Jx2 

t 



Yi,^i .^.^ (is + F" + G" (39) 
with 

fS — En „ fS-Sn „ 

Y'/^" = is-X2)--^dB{x2)X'^'"''' = {s-X2)--^dZ^^\s), sG[0,T](40) 

Jo Jo 

As in the proof of ([36|) and (p8] we can show that 

P I sup F" > cal ) < cxD and P I sup G" > cq^ ) < 00 
V*e[o,T] J „ \te[o,T] J 

The approximation result to X"^'^ is stated in the next proposition. We will use the process 
Y3,H,n ijjg^ga,d of Y '^^^'^ because clearly they are very close and one can replace the other. 



Proposition 5 For < max (^=r|^, < < i fixed, let Y^'"''', y^,^^," and X"^^^ 

he given by [T3\) . (34\ ) and ( f3P|) respectively. Define 



^t' 



I dsy,i'^'"y/'^'", tG[o,r]. (4i) 
^0 

Then for every 7 such that < 7 < /3 and 7 + /3 < \, 

P (li^{|lx2.^." - X^'^ll^,^ > Cn-^^l^-P~^\\ognf/^}) = 0. 
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Proof: The proof follows from the proofs of Proposition [2] and Proposition U] since for every 
s we have 

■ 

Let us summarize the conclusions of Proposition [2l [J] and [5] in the main result of our 

paper. 

Theorem 2 Let be the Rosenblatt process ^ and < max ^ 3!^^ ' 2+2H ^ ^ ^ ^ 
fixed. Define 

= Xl'^'"" + 2Xi2'^'" + Xf'"'"", t G [0, T] 

with X^'^''^,X'^'^'",X^'^''^ given by l[2S\) . (4^ , [3^) respectively. Then for every 7 such 
that < 7 < /3 and -/ + (3 < ^, 

P (lkK{||X^'" - X^lloo.T > Cn-(i/2-/5-^)(logn)5/2}) = 0. 

Remark 5 The slowest rate of convergence is obtained for H close to one because in this 
case (3 is close to ^. When H is close to ^ then (3 is close to |. But this situation cannot be 
compared with previous results in the literature because the Rosenblatt process is not defined 
forH=\. 
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